Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN I^TeX style file v2.2) 



Deuterated hydrogen chemistry: Partition functions, 
EquiHbrium constants and Transition intensities for the Hg^ 
system 
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ABSTRACT 



Hg and the deuterated isotopomers are thought to play an important role in interstel- 
lar chemistry. The partition functions of H;^ D2H+ and D^ are calculated to a tem- 
perature of 800 K by explicitly summing the ah initio determined rotation-vibration 
energy levels of the respective species. These partition functions are used to calcu- 
late the equilibrium constants for nine important reactions in the interstellar medium 
involving Hg and its deuterated isotopomers. These equilibrium constants are com- 
pared to previously determined experimental and theoretical values. The Einstein A 
coefficients for the strongest dipole transitions are also calculated. 

Key words: ISM: evolution, ISM: general, ISM: kinematics and dynamics, ISM: 
molecules, astrochemistry, molecular data. 



1 INTRODUCTION 

Deuterium chemistry in the interstellar medium has had 
renewed interest of late, this is due in part to recent ob- 
servation s of multiply deuterated species in the interstellar 
medium llLoinard ct al.""200lt rCeccarellil I200I IVastel et alJ 



120041) . The cosmic abundance of deuterium with respect to 
hydrogen at low temperat ure is approximately 1.4x10^'' in 
the solar neighbourhood, JLemoine et al.l^l999^ . but a much 
higher ratio is observed between molecules and their deu- 
terium baring isotopomers in some environments. Molecules 
containing deuterium can become highly fractionated in the 
gas phase at low temperatures. This effect is thought to be 
primarily through reactions with H2D"^. H2D^was first de- 
tected by Stark ct al. (1999) in the interstellar medium and 
more recently by ICaselli et al.l (120031). However modelling 
of interstellar deuterium chemistry bv lRoberts et alj (120031) 
suggest that that all the deuterated Hg" isotopomers have 
an effect on fractionation. That is the inclusion of Dg" and 
D2H^ enhances fractionation significantly. 

The other mechanism by which deuterium fractionation 
can occur is via active grain chemistry. That is to say that 
molecules which are stuck onto grains are de uterated by deu- 
terium which accrete onto the grain mantle. iBacmann et alJ 
l|2003l) suggest that in prestellar cores, deuteration increases 
with increasing CO depletion. This would indicate that 
deuteration is primarily the product of gas phase chemistry 
as CO depletion leads to an inc rease in [H2D^]/[H;^] abun- 
dances JDalgarno fc Lepdll984^ . 

Hydrogenic gas phase reactions involving Hjj" and the 



deuterated isotopomers that are regarded to be significant 
in gas phase deuteration are tabulated in table Q Within 
the thermodynamic equilibrium regime, the equilibrium con- 
stants, K, of these reactions can be calculated from the par- 
tition functions of the reactant and product species. It is 
not possibly to separate the forward and backward rates in 
these calculations thus comparisons can only been made to 
the ratio of the forward k/, and backward rates k^. 



kb 



(1) 



Experiments measuring both the forward and back- 
wa rd rates for the reaction s of interest have b een conducted 
by I Adams fc SmithI lHJsih.lGiles et alJ Jl992l) and most re- 
cently bv lCerlich et al.l (120021) . Both Adams and Smith and 
Giles et al. u sed a variable tem perature selected ion fiow 
arrangement iSmith et al.lll982h . while Gerlich et al. used a 
low temperature multipole ion trap. 

Previously theoretical studies have concentrated on the 
H;^ partition functions and associated equilibrium constants 
for r eactions of interest; most notably in Sidhu ct al. ( 1993) 
and iNeale fc Tennvsonl (|l99a). Sidhu et al. calculated the 
partition function of Hj and H2D"'" to a temperature of 
2800 K. Sidhu et al. also calculated a number of equilib- 
rium constants. Neale et al. calculated the high tempera- 
ture H^ partition function which extended to lOOOOK. We 
recalculate the Repartition function in order to compare 
with these pre vious calculations and veriiy our methodology. 
iHerbsti Jl983) calculated partition functions using what was 
then rather limited experimental and calculated molecular 
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data, from these equilibrium constants were determined in 
a similar manner to this work. 

Many attempts have been made to observe the deuter- 
ated species of H;^. To this end the Einstein A coefficients 
of the dipole transitions have been calculated. The strongest 
transitions have been selected and are tabulated here. 



2 METHOD 

2.1 Partition function 

The partition function of H2D+ was taken from Sidhu et al. 
The internal partition functions, Zi„t for Hj^, D;^ and D2H''" 
were computed by explicitly summing the series: 

Zint — 



^{2J + l)giexp(-- 



T 



(2) 



Table 1. Calculated enthalpies for the reactions of interest 



Reaction 



Enthalpy / K 



(a) 


H++ D -^ H2D++ H 


-597.8 


(b) 


H++ HD -^ H2D++ Ha 


-231.8 


(c) 


H2D++ HD -^ D2H++ H2 


-187.2 


(d) 


D2H++ HD -^ D++ H2 


-233.8 


(e) 


H++ D2^ H2D++ HD 


-153.0 


(f) 


H++ D2^ D2H++ H2 


-340.2 


(g) 


H2D++ D2^ D2H++ HD 


-108.4 


(h) 


H2D++ D2-> D++ H2 


-342.2 


(i) 


D2H++ D2^ D++ HD 


-155.0 



Table 2. Constants in cm ^ used to calculate d iatomic partition 
functions, taken from iHuber fc Herzberg 11979^ 



where J is the rotational quantum number, gt is the nuclear 
spin degeneracy factor for state i, C2 is the second radiation 
constant and Ei is the associated energy level relative to 
the J=0 ground state in cm~^. No distinction was made 
between rotational and vibrational energy levels. Hg" has 
only one bound electronic bound state, thus there was no 
electronic contribution to the partition functions. All the 
H^ energy levels were taken relative to the J=0 vibrational 
ground state. 

The full partition function, ztot can be written as 



Ztot 



trans 
Z Zint 



where Zi„t is the internal partition function, z*' 
translational contribution to the partition function which 
can be estimated using the perfect gas model. As all the 
reactions considered in this work conserve the number of 
particles in the system, the ratio of their translati onal par- 
tition functions is given by a simple mass factor JHavmarJ 
I1967D 



trans trans 
-,trans ^trans 



/memo 
KmArriB 



3/2 



(4) 



where mx is the mass of species X. 

2.2 Equilibrium constants 

For reaction 
A+B^C+D 



(5) 



the temperature dependent equilibrium constant, K(T), was 
calculate using the following. 



^C D 

K = , „ exp I 



kTj 



(6) 



where ztot is the partition function incorporating transla- 
tional motion and U is the enthalpy of the reaction. The 
enthalpy of the reaction was calculated using 



U = Ek 



So" 



E^ 



E, 



(7) 



where E^ is the zero point energy of species X as measured 
on an absolute energy scale. Thus the enthalpies for reac- 
tions (b) to (i) were calculated in this way. The diatomic 
zero point energies were calculated using the constants of 
iHuber fc Herzberd il979l) (table ^ and equation [TTl The 



H2 



HD 



B, 


60.853 


30.443 


45.655 


Cte 


3.06 


1.0786 


1.986 


D, 


0.0471 


0.01141 


0.02605 


i^e 


4401.21 


3115.50 


3813.15 


UJeXe 


121.33 


61.82 


91.65 


ge 


1 


6 


6 


go 


3 


3 


6 


DS 


36118 


36406 


36749 


l^max 


14 


20 


16 



(3) 
; iu zero point energ i es of H j' and isotopomers were taken from 



iRamanlal et al.l l|200,^^ . The ground state for HjJ" is forbid- 
den by the Pauli principle; The lowest state ^ J=l, K=l, lies 
some 64.123 cm"'^ above this ground state (IRamanlal et alJ 



l2003l) . Thus for H3 the so called "rotational zero point en- 
ergy" was used, which is 4425.823 cm~^. For reaction (a) the 
difference in i onisation en ergy between H and D was taken 
to be 46.4 K |Keii3|l983). 

The diatomic partition functions needed for the equi- 
librium constants were calculated using the formulae given 
below. 



^ = ^(2J + 1) 



I J cxp 



F,^ + G,^ — Go 
fcT 



where 

F, = B,.J(J+l)-DeJ^(J+l)', 

G„ = LUe{l' + I) — i^eXeif + ^) 



(8) 

(9) 
(10) 

(11) 



The constants used were taken from iHuber fc Herzberd 
ljl97S'l and are tabulated in table |21 



2.3 Energy levels 

Energy levels for H^ , Dj and D2 H"'" were determi n ed us 



ing the DVR3D program suite of iTennvson et alJ (|2fl 
and the ultra-high accura cy ab initio procedure of 
IPolvanskv fc TennvsorJ (|l999|) which was based on the elec- 
tronic structure calculations of lCencek et alJ l|1998r . The pa- 
rameters used for th e DVR3D program suit e are t he same 
as those outlined in lPolvanskv fc TennvsorJ Jl999r .A total 
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of 40x(J+l) levels were computed for each J, up to J=14. 
This procedure gave at least 19119 rotation vibration energy 
levels for each molecule and ensured that all energy levels 
up to 10000 cm~^ were included. 

D2H+ has C2iy symmetry, this symmetry is fully rep- 
resent in the DVR3D program by the parity of the basis, 
even and odd; which means that energies with even (ge =3) 
and odd (go =6) parity are easily easily identified. H^ and 
DJ have T>3h symmetry, this symmetry is not fully repre- 
sented by the DVR3D program. Thus energies with E, Ai 
and A2 cannot be so easily identified. The Ai and A2 states 
are represented by DVR3D and have even and odd basis par- 
ity respectively. The E symmetry energies are determined by 
the fact that they are degenerate across even and odd basis 
by parities. Thus the E symmetry states can be identified 
by examining by hand the complete list of energy levels for 
both even and odd basis parities. For Hj the nuclear degen- 
eracy factors are 2, and 4 for E, Ai and A2 respectively. 
For DJ the the nuclear degeneracy factors are 8, 10 and 1 
for E, Ai and A2 respectively. 

There has been some confusion in the conventions used 
for nuclear spin degeneracy factor, g. In general astronomers 
use the nuclear spin degeneracy divided by its value for the 
dissociated atoms, while physicists use the nuclear spin de- 
generacy itself. Thus for example astronomers would take g 
for the E states of D J to be 8/27 while we have used 8. Both 
conventions are equally valid; however it is important to be 
consistent when calculating equilibrium constants that the 
same convention has been used in calculating the partition 
functions for all the species in the reaction. Throughout this 
work the "physicists" convention is used. 



Table 3. Calculated internal partition functions as a function of 
a temperature. Powers of ten given in parenthesis 



T{K) 


Hj 


D3+ 


D2H+ 


5 


5.8325{-08) 


10.0022 


6.0008 


10 


6.3491(-04) 


10.2351 


6.1281 


20 


0.0826 


12.5729 


7.8853 


30 


0.4654 


16.3763 


11.1196 


40 


1.1481 


21.1870 


15.0916 


50 


2.0197 


26.9604 


19.5932 


60 


2.9948 


33.6025 


24.5777 


70 


4.0243 


40.9934 


30.0255 


80 


5.0826 


49.0277 


35.9154 


90 


6.1579 


57.6247 


42.2237 


100 


7.2457 


66.7263 


48.9256 


150 


12.9039 


118.6694 


87.6098 


200 


19.0975 


179.7917 


133.6128 


300 


33.3860 


325.2799 


243.7791 


400 


50.0599 


498.7918 


375.2491 


500 


68.9696 


701.1151 


527.2953 


600 


90.1816 


936.5854 


701.5534 


700 


113.9581 


1211.2871 


901.0749 


800 


140.6995 


1531.8784 


1129.5545 




H3+ 


D3+ 


D2H+ 


ao 


-35.10102 


-0.388363 


-0.975341 


ai 


72.2463 


5.65495 


7.92203 


a2 


-66.3543 


-8.53925 


-13.5834 


as 


35.3938 


5.83071 


11.0576 


a4 


-11.3756 


-1.74965 


-4.45541 


as 


2.06118 


0.205985 


0.890251 


as 


-0.160957 -0.00289176 


-0.0704028 



3 RESULTS AND DISCUSSION 



3.1 Partition function 



Table 4. Fitting coefficients for the polynomial fit reo nation 1121 
of the partition functions in the temperature range 5 K to 800 K. 



Table E] present the values obtained by the explicit summa- 
tion of equation|5| It was found that at a temperature of 800 
K the inclusion of the J=14 energy levels only contributed 
only 0.02%, 0.77% and 0.35% to the internal partition func- 
tions for HJ and D;^ and D2H^ respectively. This the par- 
tition functions are valid up to a temperature of 800 K. 

Comparing the H^^ partition functions of this work and 
that of Sidhu et al. (figure we see that the two works are 
in good agreement. There is some minor disagreement at 
higher temperatures where in an y case the much more com- 
prehensive partition function of iNeale fc TennvsonI (|l993) 
should be used. 

The partition functions have been fit ted to the stan- 
dard formula, see ISauval fc TatumI il984ft . in the temper- 
ature range 5 K to 800 K using the data in table El The 



coefficients a„ 
table m 



for H^ 



D3+ 



logio(2) = X! ""(^°gio ^)"' 



and D2H"'" are tabulated in 



(12) 



Our fit is never more than 1.35%, 0.78%, and 
1.22% from the calculated values of z for Hj, D^^and 
D2H^respectively. 



3.2 Equilibrium Constants 

Table gives the logarithm equilibrium constants, K, for 
the reactions tabulated in table Q as a function of temper- 
ature. The equilibrium constants were calculated using the 
partition functions previou sly discussed . The H 2D^partition 
functions were taken from lSidhu et all Iil993)- 

There have been a few experiments where both the for- 
ward and backward reaction rates of interest have been mea- 
sured; so that the equilibrium constant may be deduced for 
comparison. A comparison of the experimental data avail- 
able to date is given in tables |S1 17| and |H| Our calculations 
generally show approximate agreement with the experimen- 
tal data. The notable exception is the recent experiment of 



iGerlich et al. (2002) for the H++ HD^H2D+-HH2, which 
disagrees with our calculations by 12 orders of magnitude. 
Gerlich et al. measured the forward and backward rates at 
a low temperature, 10 K, using a low temperature multi- 
pole ion trap. There have been no other experiments carried 
out at this low temperature, thus no direct experimental 
comparison can be made. Smith and Adams using standard 
extrapolation give an equilibrium constant of 1.7x10"*"® at 
10 K, which is in better agreement to our own result. Most 
models use the the equilibrium constants of lAdams fc SmithI 
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N _2 



O 




Figure 1. Hg partition function, z a s a function of tem perature, T. Crosses, t his calculation; dashed cu rve, fit of equation 1121 to our 
calculated data; Circles, calculation of lSidhu et alj ll993); Pluses, calculation of lNeale &: TennvsonI ll99a). 



ljl98lf) . The previously calculated value of Sidhu et al. used 
an enthalpy of 139.5 K, this enthalpy does not take into ac- 
count the previously discussed rotational zero point energy. 
Thus if an enthalpy of 231.8 K is used then an equilibrium 
constant of T.lxlO"''^'^, which is more consistent to our own 
is obtained. 

A comparison o f our equilibrium constants with those 
of Idles et al.l lll992ll for the reactions of interest are shown 
in table |S| Giles et al. give relative errors of ± 15% for the 
equilibrium constants. This relates to a much lower absolute 
value on the measured rate; Giles et al. state that certain 
systematic errors will cancel when taking the ratio of rates, 
thus producing lower error bounds. This may be rather op- 
timistic. In general there is better agreement between our 
calculation and the experiment of Giles et al. at the higher 
temperature of 300 K. This is most likely due to the less 
demanding nature of measuring the reaction rates as higher 
temperature. Giles et al. also calculate the equilibrium con- 
stants by calculating the partition functions of the reactant 
and product species, and also the enthalpy change for the 
reaction. These partition functions were obtained by explic- 



itly summing the energy levels (as equation IJJ as given by 
using the rigid rotor approximation and the relevant exper- 
imentally determined rotational constants. The rigid rotor 
approximation is problematic for the HJ system and defi- 
nitely inferior to our own ab initio energy level calculations. 
However we are generally in better agreement with Giles et 
al. theoretical equilibrium constants than their experimen- 
tally derived equilibrium constants. 

A comparison of equilibrium constants for the reaction 
Ht-|- HP -^ H2D"'"-|- H2 for a number of temperatures to 



those of lAdams fc SmithI £1^ and lHerbsJ I^19^j:3l are given 
in table |7| Adams and Smith estimate their errors on the 
reaction rates to be ±20%, this gives an error on the equi- 
librium constants of ±30%. W e are generally in good agree- 
ment with Adams and Smith. iHerbsli il982f) calculated the 
reaction constants using his calculated partition functions 
and reaction enthalpy. The partition functions for H^ and 
H2D^are determined from explicitly doing the sum as in 
equation|21 The energy levels a re fou nd b y using the spec tro- 
scopic constants of Oka (19 8(l.ll98ll) andlCarneyl lll98(t re- 
spectively. The constants of JHuber fc Herzberd il979f) were 
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Log,nT 



Figure 2. D2H"'" and Dg partition functions, z, as a function of temperature, T. Diamonds, D2H+ calculation; dashed curve, fit of 
eauation ll2l to D2H+ data; Crosses, Dg calculation; solid curve, fit of eauation ll2l to Dg data. 



used to calculate the partition functions of HD and H2. Our 
calcul ations are in good agreement with those of iHerbsd 
J1982D. 



4 CONCLUSIONS 



3.3 Transition Intensities 

Using wavefunctions produced in the calculation of the en- 
ergy levels outlined in section 12.31 and additional wave- 
functions calculated for H2D^ in the sam e manner; the 
DIPOLE3 program of)Tennvson et al.lll2004l : and the dipole 
surface of iRohse et alJ (|l99^; the transition intensities 
were calculated. The dipole transition intensities for H2D"^, 
D2H+and D^are given in tables lUl IIUI and 1111 respectivelv: 
transitions are given up to J=5 and a maximum frequency of 
5000 cm"'^; transitions whose relative intensity is less than 
0.0001 are neglected. The Dj energy levels are labelled by 
the notation {ui, U2, J, G, U) iWatson et al. (198;^ . The 
quantum numbers v\, vi, G and U wer e assigned by refer- 
ring to the work of lAmano et al.l il994f) and by inspection. 
The H2D"'"and D2H"'"levels are assigned by hand using the 
standard quantum numbers J, Ka, and Kc. 



Partition functions for Hg D2H 



and D3 have been cal- 
culated. These have been used to calculate the equilibrium 
constants for the important gas phase reactions involving 
Hg" and the deuterated isotopomers. In general ah initio 
calculations are easier to perform at very low temperatures 
because they require few levels whereas these very low tem- 
perature experiments are extremely challenging. In addition 
transition data have been calculated which should aid in the 
observation of deuterated isotopomers of H J . It is hoped our 
data facilitate the understanding of deuterium chemistry in 
the interstellar medium. 
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Table 5. Equilibrium constants, K, for tiie reactions given in table^ as a function of temperature. The first line corresponds to reactants while the seco 
corresponds to products. 



T (K) H++D H++HD H2D++HD 

H2D++H H2D++H2 D2H++H2 



D2H++HD 
D++H2 



H2D++HD 



H++D2 
D2H++H2 



H2D++D2 
D2H++HD 



H2D++D2 
D++H2 



D2H++D2 
D++HD 



5 


59.1951 


26.5161 


16.1410 


19.6059 


20.5256 


36.6666 


10.1505 


29.7564 


13.6154 


10 


29.2007 


12.4148 


8.0189 


9.4544 


9.8455 


17.8644 


5.4496 


14.9040 


6.8851 


20 


14.1069 


5.2661 


4.0622 


4.3690 


4.4023 


8.4645 


3.1984 


7.5674 


3.5052 


30 


9.3591 


3.1640 


2.5237 


2.6552 


2.8620 


5.3857 


2.2218 


4.8770 


2.3533 


40 


7.0102 


2.1462 


1.7791 


1.8049 


2.1178 


3.8969 


1.7507 


3.5556 


1.7765 


50 


5.6306 


1.5816 


1.3426 


1.3131 


1.7047 


3.0473 


1.4657 


2.7787 


1.4361 


60 


4.7257 


1.2337 


1.0659 


0.9998 


1.4456 


2.5115 


1.2777 


2.2776 


1.2117 


70 


4.0858 


0.9997 


0.8802 


0.7841 


1.2668 


2.1471 


1.1473 


1.9314 


1.0512 


80 


3.6116 


0.8337 


0.7465 


0.6256 


1.1377 


1.8842 


1.0505 


1.6760 


0.9296 


90 


3.2462 


0.7093 


0.6456 


0.5030 


1.0398 


1.6854 


0.9761 


1.4792 


0.8335 


100 


2.9563 


0.6122 


0.5664 


0.4045 


0.9632 


1.5296 


0.9174 


1.3219 


0.7555 


150 


2.1009 


0.3270 


0.3283 


0.0973 


0.7441 


1.0724 


0.7454 


0.8426 


0.5143 


200 


1.6813 


0.1816 


0.2021 


-0.0694 


0.6415 


0.8436 


0.6620 


0.5925 


0.3904 


300 


1.2660 


0.0311 


0.0676 


-0.2472 


0.5457 


0.6133 


0.5822 


0.3350 


0.2674 


400 


1.0602 


-0.0445 


-0.0014 


-0.3386 


0.5062 


0.5048 


0.5493 


0.2108 


0.2122 


500 


0.9374 


-0.0893 


-0.0421 


-0.3927 


0.4914 


0.4493 


0.5387 


0.1460 


0.1880 


600 


0.8569 


-0.1179 


-0.0687 


-0.4274 


0.4904 


0.4217 


0.5397 


0.1123 


0.1810 


700 


0.7999 


-0.1380 


-0.0866 


-0.4507 


0.4967 


0.4101 


0.5481 


0.0974 


0.1840 


800 


0.7579 


-0.1523 


-0.0997 


-0.4668 


0.5076 


0.4079 


0.5602 


0.0934 


0.1931 



Table 6. A comparison of Equilibrium constants with lGiles et alj il99 



80K 



300K 





Giles et al. 


Giles et al. 


This work 


Giles et al. 


Giles et al. 


This work 




(Expt) 


(Theory) 




(Expt) 


(Theory) 




H+-)- HD -► H2D+-I- H2 


3.8 (±0.6) 


6.6 




6.82 


1.80 (±0.3) 


2.0 




1.07 


H2D+-I- HD -^ D2H+-I- H2 


1.7 (±0.3) 


3.6 




5.58 


0.80 (±0.1) 


0.8 




1.17 


D2H+-I- HD -^ D+-I- H2 


1.8 (±0.3) 


2.0 




4.22 


0.40 (±0.1) 


0.3 




1.00 


H+-I- D2^ H2D+-I- HD 


8.2 (±1.2) 


13.2 




13.73 


5.20 (±0.8) 


6.7 




3.51 


H+-F Da^ D2H+-H H2 


9.2 (±1.4) 


48.3 




76.59 


5.30 (±0.8) 


5.1 




4.10 


H2D++ D2^ D2H+-I- HD 


4.4 (±0.7) 


7.2 




11.23 


1.90 (±0.3) 


2.5 




3.82 


H2D+-I- D2^ D+-I- H2 


4.4 (±0.7) 


14.5 




47.43 


0.70 (±0.1) 


0.8 




2.16 


D2H+-H D2^ D+-h HD 


1.7 (±0.3) 


4.0 




8.50 


0.80 (±0.1) 


1.0 




1.85 



Table 7. A comparison of Equilibrium constants with I Adams fc Smitbl il98U) and iHerbstI (l 
H2. 



for the reaction Hg + HD 



H2D^ 



T(K) Adams and Smith Herbst 



This work 



80 4.48 (±1.3) 5.9 

200 2.35 (±0.7) 2.6 

295 1.96 (±0.6) 2.1" 



6.82 
1.52 
1.07" 



a. The theoretical value is actually at 300 K 



Table 8. A comparison of Equilibrium constants at a temperature of 10 K for the reaction Hl|'+ HD -^ H2D++ H2. Powers often giv 
in parenthesis. 



This Work 2.6(+12) 

Gerlich et al. Gerlich et al. ('2002') 7.14 

Adams and Smith Adams fc Smith C1981') 1.5(+9)" 

Sidhuetal.? 7.1(+12)'' 

a. This value is extrapolated from experimental data 

b. This value uses the corrected enthalpy of 231.8 K 
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Table 9. Einstein A coefficients for transitions form fow-lying levels of H2D"'". Powers of ten given in parenthesis. 



J' 


< 


K 




e' / cm-l 


j" 


K 


< 




e" / cm-i 


a;ij(calc.) / cm ^ 


uJij-(ohs.) 1 cm -•- 


A,/ / s-i 


1 


1 





Bi 


72.457 




1 




Ai 


60.027 


12.429 


- 


1.2186(-4) 


1 





1 


A2 


45.698 









Ai 


0.000 


45.698 


- 


4.0397(-3) 


2 


1 


2 


A2 


138.843 




1 




Ai 


60.027 


78.816 


- 


1.8762(-2) 


2 





2 


Ai 


131.638 









A2 


45.698 


85.94 


- 


3.0338(-2) 


2 


1 


1 


Bi 


175.939 




1 




Bi 


72.457 


103.483 


- 


4.238{-2) 


2 


2 





Ai 


223.868 









A2 


45.698 


178.17 


- 


1.6643(-2) 











Ai 


2205.916 









A2 


45.698 


2160.218 


2160.176° 


17.545 


1 


1 





Bi 


2278.465 




1 




Ai 


60.027 


2218.438 


2218.393" 


10.372 


1 





1 


A2 


2246.727 









Ai 


0.000 


2246.727 


2246.697'" 


1.8962 


2 





2 


Ai 


2318.377 









A2 


45.698 


2272.68 


- 


0.35157 











A2 


2335.338 




1 




Ai 


60.027 


2275.31 


2275.403'" 


145.65 


1 





1 


Ai 


2383.878 




1 




Bi 


72.457 


2311.421 


2311.512'" 


83.419 


1 


1 





Bi 


2409.227 









A2 


45.698 


2363.529 


- 


78.984 


2 


2 





Ai 


2427.119 









A2 


45.698 


2381.421 


2381.367" 


3.058 


1 


1 


1 


A2 


2402.699 









Ai 


0.000 


2402.699 


2402.795'" 


60.938 


2 





2 


A2 


2477.681 




1 




Ai 


60.027 


2417.654 


2417.734'" 


29.82 


2 


1 


2 


Ai 


2490.966 









A2 


45.698 


2445.268 


2445.348'" 


58.586 


2 


2 


1 


Bi 


2568.382 




1 




Bi 


72.457 


2495.925 


2496.014'" 


60.067 


2 


2 





A2 


2569.489 




1 




Ai 


60.027 


2509.461 


2509.541'" 


48.476 











Ai 


2992.524 









A2 


45.698 


2946.826 


2946.802'' 


53.167 


1 


1 





Bi 


3063.331 




1 




Ai 


60.027 


3003.304 


3003.276'' 


27.509 


1 





1 


A2 


3038.198 









Ai 


0.000 


3038.198 


3038.177'' 


20.353 


2 


1 


2 


A2 


3128.888 




1 




Ai 


60.027 


3068.86 


3068.845'' 


20.088 


2 





2 


Ai 


3123.324 









A2 


45.698 


3077.626 


3077.611'' 


24.757 


2 


1 


1 


Bi 


3167.147 




1 




Bi 


72.457 


3094.69 


3094.671'' 


19.302 


2 


2 





Ai 


3209.847 









A2 


45.698 


3164.149 


- 


1.5976 











Ai 


4287.61 









A2 


45.698 


4241.912 


- 


16.953 


1 





1 


A2 


4331.45 









Ai 


0.000 


4331.45 


- 


9.6306 


2 


1 


2 


A2 


4412.461 




1 




Ai 


60.027 


4352.434 


4352.360'= 


14.871 


2 





2 


Ai 


4407.925 









A2 


45.698 


4362.227 


- 


15.45 











A2 


4461.832 




1 




Ai 


60.027 


4401.805 


- 


88.628 


2 


1 


1 


A2 


4512.486 









A2 


45.698 


4466.788 


- 


0.41306 


1 


1 





Bi 


4536.348 









A2 


45.698 


4490.65 


- 


44.342 


2 





2 


A2 


4555.88 




1 




Ai 


60.027 


4495.853 


4495.881'= 


27.013 


1 


1 


1 


A2 


4512.558 









Ai 


0.000 


4512.558 


4512.567'= 


41.115 


2 


1 


2 


Ai 


4563.405 









A2 


45.698 


4517.707 


- 


38.458 











Ai 


4602.746 









A2 


45.698 


4557.048 


- 


69.044 


2 


2 


1 


Bi 


4677.548 




1 




Bi 


72.457 


4605.091 


- 


38.831 


1 


2 


1 


Bi 


4677.74 




1 




Ai 


60.027 


4617.713 


- 


33.563 


2 


2 





A2 


4691.531 




1 




Ai 


60.027 


4631.504 


- 


20.198 


1 





1 


A2 


4657.859 









Ai 


0.000 


4657.859 


- 


9.0071 


2 





2 


Ai 


4761.399 









A2 


45.698 


4715.701 


- 


7.7521 


2 


2 





Bi 


4845.211 









A2 


45.698 


4799.514 


- 


1.3349 











Ai 


5039.84 









A2 


45.698 


4994.142 


- 


10.723 


a. 


Frequencies of 
Frequencies of 
Frequencies of 


Foster et al.1 Jl986D 


















b. 


Kozin et alj ll988r 




c. 


Farnik et al.l J2002l 





REFERENCES 

Adams N. G., Smith D., 1981, Astrophysical Journal, 248, 

373 
Amano T., Chan M., Civis S., McKellar A., Majewski W., 

Sadovsii D., Watson J., 1994, Can J. Phys., 72, 1007 
Bacmann A., Lefloch B., Ceccarelli C, Steinacker J., 

Ccistets A., Loinard L., 2003, Astrophysical Journal, 585, 

L55 
Carney G. D., 1980, Chem. Phys., 54 
Caselli P., van der Tak F. F. S., Ceccarelli C., Bacmann 

A., 2003, Astro. Astrophys., 403, L37 
CeccareUi C, 2002, Planet Space Sci., 50, 1267 



Cencek W., Rychlewski J., Jaquet R., Kutzelnigg W., 1998, 
J. Chem. Phys., 108, 2831 

Dalgarno A., Lepp S., 1984, Astrophysical Journal, 287, 
L47 

Farni'k M., Davis S., Kostin M. A., Polyansky O. L., Ten- 
nyson J., Nesbitt D. J., 2002, J. Chem. Phys., 116, 6146 

Foster S. C, McKellar A. R. W., Peterkin I. R., Watson 
J. K. C, Pan F. S., 1986, J. Chem. Phys., 84, 91 

Gerlich D., Herbst E., Roueff E., 2002, Planetary and Space 
Science, 50, 1275 

Giles K., Adams N. G., Smith D., 1992, Journal of Physical 
Chemistry, 96, 7645 



8 Jayesh Ramanlal and Jonathan Tennyson 

Table 10. Einstein A coefficients for transitions form low-lying levels of D2H+. Powers of ten given in parenthesis. 



J' 


< 


< 




e' / cm-i 


j" 


< 


< 


] 


e" / cm-i 


u)if{ca\c.) 1 cm ^ 


ti.'i/(obs.) / cm -"- 


A,^ / s-i 


1 


1 





Bi 


57.993 









Ai 


34.918 


23.075 


- 


5.0911(-4) 


1 


1 


1 


A2 


49.255 









Ai 


0.000 


49.255 


- 


3.3031(-3) 


2 





2 


Ai 


101.716 




1 




A2 


49.255 


52.461 


- 


2.0702(-3) 


2 


1 


2 


A2 


110.259 









Ai 


34.918 


75.341 


- 


1.0597(-2) 


2 


2 


1 


Bi 


179.173 




1 




Bi 


57.993 


121.18 


- 


4.4504(-2) 


2 


2 





Ai 


182.074 




1 




A2 


49.255 


132.819 


- 


4.4629(-2) 











Ai 


1968.146 




1 




A2 


49.255 


1918.89 


1918.908° 


52.175 


1 





1 


Ai 


1998.523 




1 




Bi 


57.993 


1940.53 


1940.551° 


21.989 


1 


1 





Bi 


2027.034 









Ai 


34.918 


1992.116 


1992.130° 


28.48 


2 





2 


Ai 


2055.077 




1 




A2 


49.255 


2005.821 


2005.844° 


6.0891 


1 


1 


1 


A2 


2014.09 









Ai 


0.000 


2014.09 


2014.106° 


14.364 


2 


1 


2 


A2 


2062.923 









Ai 


34.918 


2028.005 


2028.024° 


9.9809 











A2 


2078.435 









Ai 


34.918 


2043.517 


2043.515° 


21.728 


1 


1 


1 


Ai 


2128.7 




1 




Bi 


57.993 


2070.707 


2070.708° 


16.782 


1 


1 





A2 


2136.248 




1 




A2 


49.255 


2086.992 


2086.990° 


12.135 


2 


2 


1 


Bi 


2145.612 




1 




Bi 


57.993 


2087.619 


2087.630° 


11.233 


2 


2 





Ai 


2149.555 




1 




A2 


49.255 


2100.299 


2100.307° 


5.9919 


1 





1 


A2 


2118.589 









Ai 


0.000 


2118.589 


2118.588° 


14.303 


2 


1 


2 


Ai 


2202.779 




1 




A2 


49.255 


2153.524 


2153.525° 


21.901 


2 





2 


A2 


2194.064 









Ai 


34.918 


2159.146 


2159.145° 


19.827 


2 


1 


1 


Bi 


2225.161 




1 




Bi 


57.993 


2167.168 


2167.166° 


17.293 


2 


2 





A2 


2257.594 









Ai 


34.918 


2222.675 


- 


0.38661 











Ai 


2736.98 




1 




A2 


49.255 


2687.724 


- 


86.183 


1 





1 


Ai 


2771.523 




1 




Bi 


57.993 


2713.53 


- 


44.856 


1 


1 





Bi 


2793.96 









Ai 


34.918 


2759.042 


2759.036'' 


44.572 


1 


1 


1 


A2 


2785.338 









Ai 


0.000 


2785.338 


2785.332*' 


31.182 


2 





2 


Ai 


2837.556 




1 




A2 


49.255 


2788.3 


2788.300'' 


16.939 


2 


1 


2 


A2 


2845.72 









Ai 


34.918 


2810.802 


2810.800'' 


29.412 


2 


2 


1 


Bi 


2912.708 




1 




Bi 


57.993 


2854.716 


2854.707'' 


29.321 


2 


2 





Ai 


2915.616 




1 




A2 


49.255 


2866.36 


2866.350'' 


22.136 











Ai 


3821.309 




1 




A2 


49.255 


3772.054 


- 


6.5942 


1 





1 


Ai 


3851.977 




1 




Bi 


57.993 


3793.984 


- 


5.0848 


1 


1 





Bi 


3881.727 









Ai 


34.918 


3846.809 


3846.786'= 


3.3258 


2 





2 


Ai 


3909.933 




1 




A2 


49.255 


3860.677 


3860.660"= 


3.0475 


1 


1 


1 


A2 


3871.398 









Ai 


0.000 


3871.398 


3871.377"= 


3.6663 


2 


1 


2 


A2 


3921.988 









Ai 


34.918 


3887.07 


3887.052"= 


5.4261 


2 


2 


1 


Ai 


4010.528 




1 




Bi 


57.993 


3952.535 


- 


3.6408 


2 


2 





Ai 


4013.18 




1 




A2 


49.255 


3963.924 


- 


3.4514 











Ai 


4042.815 




1 




A2 


49.255 


3993.56 


3993.518"= 


47.457 


1 





1 


Ai 


4058.521 




1 




Bi 


57.993 


4000.528 


4000.494"= 


42.9 











A2 


4060.822 









Ai 


34.918 


4025.904 


4025.873"= 


41.465 


2 





2 


Ai 


4097.094 




1 




A2 


49.255 


4047.839 


4047.840"= 


21.275 


1 


1 


1 


A2 


4062.925 









Ai 


0.000 


4062.925 


4062.889"= 


29.138 


2 


1 


2 


Ai 


4097.933 









Ai 


34.918 


4063.015 


4062.983"= 


26.675 


1 


1 





A2 


4101.122 









Ai 


34.918 


4066.204 


4066.158"= 


23.775 


1 





1 


Bi 


4119.147 




1 




A2 


49.255 


4069.891 


4069.859"= 


20.300 


2 


2 


1 


Bi 


4179.804 




1 




Bi 


57.993 


4121.811 


- 


24.474 


1 


1 


1 


A2 


4122.993 









Ai 


0.000 


4122.993 


- 


1.3689(-3) 


2 





2 


Ai 


4214.033 




1 




A2 


49.255 


4164.VVV 


- 


5.7652 


2 


2 





A2 


4208.006 









Ai 


34.918 


4173.088 


- 


0.92222 


2 


1 


2 


Ai 


4229.853 




1 




A2 


49.255 


4180.597 


- 


2.4239 


2 


2 


1 


A2 


4252.4 









Ai 


34.918 


4217.482 


- 


1.3913 











Ai 


4648.808 




1 




A2 


49.255 


4599.553 


- 


10.723 


1 





1 


Ai 


4673.469 




1 




Bi 


57.993 


4615.476 


- 


9.1384 











A2 


4674.96 









Ai 


34.918 


4640.042 


- 


13.692 


2 





2 


Ai 


4720.562 




1 




A2 


49.255 


4671.307 


- 


4.6985 


1 


1 





Bi 


4706.784 









Ai 


34.918 


4671.866 


- 


4.9443 


1 


1 


1 


A2 


4681.85 









Ai 


0.000 


4681.85 


- 


7.1377 


1 


1 





Ai 


4732.173 




1 




A2 


49.255 


4682.918 


- 


6.9989 


2 


1 


2 


A2 


4723.648 









Ai 


34.918 


4688.73 


- 


6.6333 


1 





1 


A2 


4727.065 









Ai 


0.000 


4727.065 


- 


0.44733 


2 


2 


1 


Bi 


4796.924 




1 




Bi 


57.993 


4738.931 


- 


6.7834 


2 


2 





Ai 


4798.798 




1 




A2 


49.255 


4749.542 


- 


2.1769 


2 





2 


A2 


4807.782 









Ai 


34.918 


4772.864 


- 


0.87466 


2 


n 


n 


Ri 


4852.274 









Ai 


34.918 


4817..35fi 


_ 


0.5031 1 



Deuterated hydrogen chemistry: Partition functions, Equilibrium constants and Transition intensities for the H^ s\ 



Table 1 1 . Einstein A coefficients for transitions form low- lying levels of DJ 



G U 



E_/ 



G 



E / LOi j{ca[c.) I tJij(obs.)" / 









01 


1 




E 


1834.655 




00 







e 


32.324 


1802.331 


1802.349 


264.18 


1 


01 





-1 


a;' 


1884.380 




00 







Al' 


43.609 


1840.770 


1840.789 


34.99 


1 


01 


1 




E' 


1878.561 




00 







e" 


32.324 


1846.237 


1846.256 


141.96 


1 


01 







A'; 


1888.048 




00 







a; 


0.000 


1888.048 


1888.065 


252.84 


2 


01 


1 


-1 


e' 


1955.981 




00 







e" 


32.324 


1923.657 


1923.670 


107.50 


2 


01 







a; 


1979.203 




00 







i* 


43.609 


1935.593 


1935.609 


24.65 


2 


01 


1 




E 


1968.077 




00 







32.324 


1935.753 


- 


119.38 


1 


02 


3 


2 


a; 


3646.285 




00 







E 


43.609 


3602.676 


3602.669 


234.05 





02 


2 


2 


E 


3650.700 




00 







32.324 


3618.376 


3618.371 


231.47 


2 


02 


4 


2 


e' 


3662.342 




00 







E" 


32.324 


3630.018 


3630.022 


140.29 


1 


02 


-3 


-2 


a;' 


3647.190 




00 







a; 


0.000 


3647.190 


- 


193.96 


1 


02 


2 


2 


e' 


3694.893 




00 







e" 


32.324 


3662.569 


3662.557 


114.64 


2 


02 


-3 


2 


a;: 


3736.400 




00 








a; 


43.609 


3692.791 


3692.785 


17.29 


2 


02 


2 


2 


e' 


3783.234 




00 







e" 


32.324 


3750.909 


3750.879 


23.66 



a. Frequencies of Amano et al. lAmano et al.l il994^ 
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